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I was asked by the organizers of "Strings to LHC" to provide a brief overview of brane cos- 
mology. This is a large subject, and my scope will be necessarily limited to some of the key 
ideas. In particular, I will not discuss the many interesting developments with higher-codimension 
branes, rather focusing on the simplest case of codimension-one. Furthermore, this is not a review 
of string cosmology, nor of brane-antibrane inflation. Braneworld cosmology is a string-inspired 
framework, where the "branes" are merely delta-function sources in the stress-energy tensor, which 
generally do not have the other properties of D-branes, such as charge and internal dynamics. I will 
give a historically ordered account stalling with large extra dimensions, then turning to warped 
geometries and the DGP model. 

1. Flat extra dimensions 

The modern era of brane cosmology began with ref. ^ in the context of large extra dimen- 
sions, made possible by the hypothesis that the standard model of particle physics is localized on a 
D-brane. Since only gravity was presumed to propagate in the extra dimensions, only short-distance 
(sub-millimeter) probes of the gravitational force could constrain the size of the extra dimensions. 
Experimentally it was known that gravity did not deviate from its known four-dimensional form 
down to distances just below a millimeter, hence allowing for extra dimensions of that magnitude. 
Additionally, with two extra dimensions, the fundamental gravity scale would be lowered to a 
value consistent with the TeV scale, using the relation 



in d spacetime dimensions, with d — 4 having size R. Cosmology becomes highly constrained in 
such a scenario because the phase space for Kaluza- Klein gravitons g is so large that they would be 
overproduced if the temperature of the universe was ever greater than some maximum value ~ 1 
MeV. The gravitons are populated by interactions like 77 — > gg. Even though each interaction 
vertex is suppressed by 1 /M p in the 4D effective theory, the phase space is so dense that the total 
cross section is similar to having couplings that are only suppressed by l/Mg, where at the 
1 — 10 TeV scale. The KK gravitons are constrained because they can overdose the universe, or 
distort the diffuse 7 ray spectrum by their decays at late times. 

One of the key papers motivating interest in brane cosmology was ref. % (BDL). Their 
work was inspired by string theory, the Hofava-Witten model |^], in which the 11th dimension 
of heterotic M-theory is compactified on an interval y G [0, 1], which can be considered as a circle 
parametrized by y G [—1,1] with points identified under a Z2 orbifold symmetry y — > —y. The ends 
of the interval are the orbifold fixed points. BDL considered a 5D line element of the form 



with branes at y = and y = 1 , and bo assumed to be constant. 

The branes provide delta-function sources in the Einstein equations, which take the simple 
form 



Mp ~ M d d - 2 R l, - A 



(1.1) 



ds 2 = -n 2 (y)dt 2 +a 2 (t,y)dx 2 +bf } dy 2 



(1.2) 




(1.3) 
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For constant equation of state /?/p = won both branes, the solution is also simple, 

boPo 
1 



l-J^\y\\ (1.7) 



2M\ V 3 



1 + S(w + ^)M (1.8) 



&0P0 

where q~ l = 3(1 + vv). The time-dependence of the scale factor differs from that in standard 
cosmology, where q~ l = 3(1 + w)/2. A simple way of seeing why this happens is to consider 
the Einstein equation (|L3|). By integrating it over the region y £ [— e, e], we find that 



Ad d bop 



_^ = _ =2 -=--^ (1.9) 
e a a a 3M| 



where we used d(—e) = —d{e) due to the Z2 orbifold symmetry. Similarly, integrating (L4) in 
this region gives 

- = — %(3p-p) (1.10) 



Using these values in ( 1.5) and choosing n = 1 at y = gives 



H2+ ; = 34 P(P " 3 '' ) <U1) 
Because of energy conservation, p = 3H(p + p), which follows from (1.6), eq. ( |1 . 1 1| ) is consistent 
with the first-order equation 

H 2 =(-^] (1.12) 

This modified form for the Friedmann equation can be solved for a(t) given the usual scaling 
p ~ a~ 3 ( 1+w \ to show that a(t) ~ t q as in (O). 



Needless to say, the dependence H ~ p instead of the standard relation // = y/SnGp /3 is 
ruled out by the successful prediction of big bang nucleosynthesis (BBN) for the light element 
abundances. BBN is quite sensitive to changes in the expansion rate of the universe. 

2. Warped extra dimensions 

To overcome the problem with BDL's Friedmann equation, a simple modification was pro- 
posed in [Q] and [Bp: add a negative bulk cosmological constant, A^, which further modifies the 
Friedmann equation to the form 

H 2 =(-P 1 )\^ (2.1) 
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This greatly improves the situation since we can let 

p = x + p m (2.2) 

where X represents the tension of the brane at y = 0, while p m stands for ordinary matter or radiation 
on the brane. By tuning X — > To in such a way that the constant contribution to H vanishes, 

T ° V + ^=0 (2.3) 



6M| J 6M] 

which is the usual tuning of the cosmological constant to zero, we are left with the modified Fried- 
mann equation 

l 

where we should identify the coefficient as 



# = rPm M + P- (2.4) 



T ° ?G (2.5) 



18Mf 3 

This is interesting because it reproduces the standard 4D result at low energies p m <C To, but predicts 
the BDL-like modification H ~ p in the high-energy regime, p m » To- 

Moreover, the background solution in the presence of < becomes the warped solution of 
Randall and Sundrum (RS) Q in the limit p,„ — > 0, 

a(y)=n{y)=e- bk \y\ (2.6) 

thus providing a link between brane cosmology and the warped compactification scenario which 
was proposed for completely noncosmological reasons — namely to solve the hierarchy problem. 
(If T is not tuned to the static value To, one finds warped de Sitter solutions a(t,y) = ^'(coshMlyl — 
(x/x )smhbk\y\) [0].) 

However this was not a complete solution to the problems which plagued the BDL model. 
In both BDL and in RS it was necessary to have a negative tension brane at y = 1, with equal 
and opposite tension to the y = brane. This is not fatal since within string theory there exist 
physically consistent negative-tension objects, orientifold planes. The real problem was that the 
matter densities on the two branes had to also be tuned, such that 

Pi = -e 2kbo Po (2.7) 

where the subscript indicates the position of the brane. String theory does not provide any sources 
of matter with negative energy density. Furthermore, it was not clear why it should be necessary to 
tune the matter densities on the two branes. 

In retrospect, it should have been obvious that something else was seriously lacking in the 
model of [0] since it failed to reproduce the predictions of general relativity (GR) even if the 



extra dimension became arbitrarily small — notice that the strange Friedmann equation (1.12) is 
completely insensitive to the value of bo, which is the size of the extra dimension. The problem 
was shown to be due to artificially assuming that b was static, even though there was no mechanism 
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in the theory to stabilize it at any particular value This criticism applied equally to the warped 
models of [|,|]. 

In fact, bo should be treated as a massless scalar field b(x^), the radion, which leads to a 
scalar-tensor theory of gravity rather than GR. Although it is possible to find solutions where this 



massless scalar is constant by tuning the brane matter densities as in ( |2.7[ ), fifth force constraints 
oblige us to give the radion a mass, and stabilize it at some particular value. Ref. [Js[] showed that 
once this is done, the standard Friedmann equation is recovered in either the warped or unwarped 
models, up to corrections which are suppressed by the radion mass, rather than the brane tension. 



It is interesting to see how the argument which led to eq. ( |1.12| ) is changed once we add a 
potential energy V(b) for the radion (and possibly a bulk cosmological constant The 5D 
Einstein equations become 

3^-3^111+1- 



ig vH£> 2 ) = w 5 (v 5{l " iy - yd)+At+v,M ) <28) 

, a n 2 f a" n" a' ( a' n'\) n 2 0/ , , A 

H 2 + 2---j 2- + - + - - + 2- H = --3 \Y,PiS(b Q (y-yi ))+A b +V(b ) \ (2.9) 
a bfi { a n a \a n ) J Mj / 



a \ n" a I a n \ n 2 



a J bn a \ a n J M 5 



— \A b + V(bo) + boV'{b )\ (2.10) 



where V(bo) is the radion potential, and we are still assuming that bo is constant in space and time 



for simplicity. The junction conditions ( |L9j]L10[ ) are unchanged by the radion potential, but the 5-5 
Einstein equation ( |2.1Q| ) no longer determines H 2 + 1 in terms of brane sources because of the extra 
V(bo) contribution. Instead this equation determines the shift in the radion due to the expansion of 
the universe, b — > bo + 8b. For the Randall-Sundrum model, ref. ^ finds 



\a 2 M 2 



where nib is the radion mass, and a is the warp factor; therefore aM p = A is the TeV (or infrared) 
scale. 

It is no longer simple to compute the modifications to the Friedmann equation once the radion 
is stabilized. Moreover the result depends upon whether the expansion of the universe is measured 
at the Planck brane (y = 0) or the TeV brane (y = 1). At the TeV brane, the leading correction was 
found to be [g] 

" 2 L, ~(p + -^(1 " 3w)(l + w) + 0(p 2 )) (2.12) 



lv=1 3 V m 2 A 2 

where c = 0(1), p is the warped energy density on the TeV brane, and w = p/p. Curiously the p 2 
correction vanishes when w = 1/3 (radiation) or w = — 1 (inflation). 



3. Randall-Sundrum II cosmology 

Despite its shortcomings (because it ignored the need for radion stabilization), the modified 
Friedmann equation (2.4) turned out to have a deeper significance. This is in connection with 
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the Randall-Sundrum II model [fTO]], in which the TeV brane is displaced to y —> oo (equivalently 
b — > oo) and thereby removed from the theory. There is no physical radion in this theory (it can 
be gauged away by a coordinate transformation), hence no need for stabilization. Eq. ( |2.4[ ) is still 
the correct Friedmann equation, but now it is fully justified, unlike in the RSI model. Big bang 
nucleosynthesis puts a weak constraint on the brane tension; r > (MeV) 4 to insure that H has the 
standard form during BBN. Terrestial measurements of the gravitational force give a much stronger 
constraint since the exchange of KK gravitons modifies the Newtonian gravitational potential to the 
form 

fflim? / 1 \ 
y= _ _^ 1 + _) (3.!, 

with k = M^/Mp. Interpreting the current limit from the Etit-Wash experiment [ |TI| ] to imply that 
kr > 1 at the scale 0.04 mm, we find the bound 



hence 



and the brane tension is 



k > 5 meV (3.2) 



M 5 = k l/3 M 2 / 3 > 3 x 10 5 TeV (3.3) 



T = 2\kM\ > (7 TeV) 4 , (3.4) 

far above the energy scale of BBN. 

3.1 Inflation with nonstandard Friedmann equation 

An obvious application of the modified Friedmann equation is to inflation at scales greater 
than t 1 / 4 , where H ~ p. Ref. [12] examined chaotic inflation in this context, noting that one can 



get inflation to work with steeper potentials than usual, since H is larger than its conventional value 
and Hubble damping is more effective. It was shown that the slow-roll parameters get modified 
relative to their conventional values (denoted by £q and rjo), 

/ i+y/T \ 4t 

£ = H a + y/2^ h y £o (3 - 5) 

1 2t 

T7 = T7 ; ~ Tin (3.6) 

' '1+V/2T V ' 

where V is the potential and the final estimates are in the limit 1^>T. One can show that inflation 
occurs when 2 < |V (p < — |p) intead of the usual condition <p < V (p < — ^p) in this regime. 
Furthermore the number of e-foldings as a function of the inflaton field is given by 

^-si/^O^) <37) 

and the scalar power amplitude is enhanced relative to its usual value, 

2 



A]=Al[l + -\ (3.8) 



v v 
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The spectral index formula keeps its usual form in the new slow-roll parameters, but their relation 
to the number of e-foldings in general differs from the standard result (shown below in parentheses) 



in the high-energy regime [ |13| ] : 



-A (c .f. _2), V = m 2 ^ 
-| (c.f. V = X^ 



-6 £ +2t7^ f v ;-;- ";z (3.9) 



Coincidentally, the prediction for 4 chaotic inflation does not change, although [|l3|] shows that in 
the intermediate region where V ~ X this is no longer the case, and the tension between X<p 4 theory 
and the CMB data is made somewhat worse. For the 2 model, the increase in \n s — 1| in the high 
V regime degrades the agreement with the data (viewed as likelihood contours in the plane of the 
tensor-to-scalar ratio r versus n s — 1) relative to the standard V <C T regime: it pushes it outside the 
1-a contour. 

A novel feature is that inflation can occur with subPlanckian field values, (p < M p , which is 
not the case for conventional chaotic inflation, and is often regarded as a drawback. However, one 
still needs (j) > M5; for V = m 2 <p 2 , inflation occurs for ~ 3OOM5, and the COBE normalization 
requires m = 5 x 10~ 5 M5. 

Braneworld inflation was also considered using the potential 

V = Voe- a */ M > (3.10) 

in [fPfl, which relies upon gravitational particle production to get reheating since there is no min- 
imum around which <p can oscillate to give ordinary reheating. In standard cosmology, one needs 
to have a 2 < 2 in order to get inflation from ( 3.10| ), but this is no longer true in the high-energy 



regime. Ref. [ 14] finds that with the COBE normalization the inflationary scale, brane tension and 



5D Planck mass are fixed in terms of a as 



10 15 GeVV /10 15 GeV\ 4 3 x 10 16 GeV 

I . t ~ I — — I . ~ 

a 



The number of e-foldings as a function of field value is given by 

N=^{V{^)-V enA ) (3.12) 

in the high-energy phase. This shows that a can be taken much larger than the standard range if T 
is appropriately adjusted downward, and so the potential can be much steeper than usual while still 
supporting slow-roll inflation. The eventual reheating temperature is found to be T ~ 10 /a GeV. 
Independently of a, the predicted relation between spectral index and number of e-foldings is 

4 

72,-1 = (3.13) 

N+l 

leading to n s = 0.92 as well as the tensor-to-scalar ratio r = 0.03. This value of n s is somewhat too 
low relative to the currently favored WMAP3 value n s = 0.95. 
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4. Randall- Sundrum I cosmology 

We now return to the more rich and complicated case of two branes in a compact warped 
geometry. Notice that eq. ( 2. 1 1[ ) predicts that the radion gets destabilized at temperature 



T ~ ^Jm b K ~ £ x TeV (4.1) 

where e ~ m^/A is a small parameter which we shall specify below. To understand the dynamics 
of this phase transition between compact and noncompact extra dimensions, we must discuss the 
mechanism of stabilization. The most popular such mechanism was proposed by Goldberger and 



Wise (GW) [15] 



4.1 Goldberger- Wise mechanism 

The GW mechanism uses a bulk scalar field <J> with Lagrangian 

Jgf = (d<S>) 2 -m 2 <S> 2 -V 8(y)-Vi8{y-\) (4.2) 

where the brane potentials Vo,i are minimized at <J> = vo and <J> = vi, respectively. There are two 
competing effects in this model. The mass term 

b J dye- 4kby m 2 d> 2 (4.3) 

is minimized when b is small, while the gradient energy 

l - J dye- 4kby & 2 (4.4) 

is minimized at large b. The brane potentials insure that the trivial solution <E> = is energetically 
disfavored; instead (assuming the coefficients of Vo,i are sufficiently large) <J> must vary between 
the boundary values vo and vi. The result of this competition is that b gets stabilized at some 
intermediate value, 

where e = m 2 /4k 2 -C 1. The canonically normalized radion field can be written approximately as 

<j>=fe- kb W, f = V6M p (4.6) 
and its VEV using ( fO| ) is related to the warp factor by 

-M=e- kh = ^l^( V A l ' E (4.7) 



/ V v o, 
Notice that the large hierarchy is naturally generated by raising a moderately small number T] to a 
moderately large power l/e. The radion potential has the approximate form 



(4.8) 



This looks like a double- well potential, fig. [T], with a slightly deeper minimum at = (0) than 
at (j> = 0. (This depth can be changed by adjusting the TeV brane tension.) The local minimum 
at (j> = corresponds to the decompactified RS II theory, where b — > oo. This effective potential 
predicts that the radion mass scales like £ 3 / 4 A [M 0], but a more careful analysis shows that 



actually my ~ eA [ 1 8 ] . 
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Figure 1: The Goldberger-Wise radion potential at zero temperature (solid curve) and at high temperature 
(dashed curve). 



4.2 Decompactification phase transition 

At high temperatures, the GW potential naturally gets a thermal correction which lifts the 
nontrivial minimum and leads to a transition between the compactified and decompactified phases 
of the theory (fig. [I]). This was first studied in Jl7|], and has been more recently considered in 
[)T9[]-[pT|]. A simple way of understanding the destabilization at high temperature is through the 
coupling of the radion to matter. It couples conformally, i.e. to the trace of the stress-energy tensor, 



This predicts no shift in the radion during a radiation-dominated era; however at high T there is 
also a correction to the radion mass, r 2 2 , and its effect is parametrically the same as that of the 
linear coupling if we simply estimate p — 3p ~ T 4 . If Vo is the zero-temperature potential, we can 
estimate the shift in the radion by taking 



Letting <p = 0o + 8<p, this implies 



d(j> 2 r A 
Using the fact that 0o = fe ~ A, we see that 

50 T 4 



aty dty A 



° 50 = m 2 h 8<$> ~ — - (4.11) 



</>o A 2 



in 



(4.12) 



b 



which becomes of order unity at the critical temperature given by eq. ( |4.1[ ). 

More insight into the high-temperature limit of the theory can be gained by considering the 
5D solution to the Einstein equations at finite temperature: this is the AdS-Schwarzschild solution 
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where there is a black hole in the bulk, located at the AdS horizon r = [22] (which corresponds 
to y = oo in the RS coordinate system), 

ds2 = -i$-^) dt2 + r2dl2+ {j2-^) dr2 (4" 13 ) 

The temperature of the system is identified with the Hawking temperature of the black hole, whose 
mass is /J., and £ = l/kis the AdS curvature length scale. In this coordinate system the branes are 
moving through the bulk, toward larger values of r, with the Planck brane at r = R(t) say, and the 
TeV brane at some smaller value of r. There is a horizon at 

r = \i l l 4 i (4.14) 

behind which the TeV brane is hidden before the compactification phase transition. 

The brane motion can be understood by considering the induced metric on the brane at r = 

R(t): 



ds=-\^--^-R 1 y-^--^J Jdt l -R z dx z = -dz 2 + R z (z)dx 2 (4.15) 

where z is the proper time for an observer on the brane. It is clear from the latter form that R(z) 
plays the role of the scale factor for the cosmological expansion measured by the brane observer. 
The dynamics of R(z) are determined by the junction condition at the brane, leading to the Fried- 
mann equation 

H 2 = ^ 4 (4.16) 

which has the interpretation of /J./R 4 being the energy density of "dark radiation," the Hawking 
radiation emitted by the black hole, consisting of thermal KK gravitons in the bulk. Notice that 
this Friedmann equation has the normal form, H 2 ~ p, rather than p 2 , since radiation redshifts 
like l/R 4 . In the AdS/CFT correspondence, these are the thermalized CFT degrees of freedom. 
This picture in which the geometry is static but the branes are moving has been called "mirage 
cosmology," [ ^3| ] since the expansion of the universe appears to be an illusion on the part of the 
brane observers, due to their motion through the bulk. 

Of course, the same physics can also be seen in the RS coordinates, where the geometry is 



explicitly time dependent and the branes are stationary [24]. The metric functions are 



a 2 (t,y) = al{t)e- 2kb y + sinh 2 ^; n(t,y) = — I ' - -^rr smh 2 kby] (4.17) 
flg(0 a V %{t) J 



and the horizon is at 



„-2kby 

e ~ Ae -Vcby_ c _ 



sinh z kby a 4 



o 

^i 1 "! < 4 - 18) 



which recedes toward the AdS horizon as the universe expands. In these coordinates c/qq is the 



dark radiation density, so c ~ }l in terms of the coordinates in (4.13). At some critical temperature 
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the TeV brane is uncloaked by the horizon, signaling the compactification phase transition [25]. In 
the prior epoch when the TeV brane is hidden below the horizon, it is as good as nonexistent, since 
it is inaccessible to observers outside the horizon. For them, the TeV brane might as well be at 
r = 0, which is the decompactified phase of the theory. 

The purely geometrical interpretation given above does not include the effects of radion sta- 
bilization. From the shape of the radion potential shown in fig. [I], we see that there is a first-order 
phase transition of a peculiar kind: within the high-J 5D single-brane universe, bubbles nucleate, 
inside of which the universe is 4D and the TeV brane has appeared. References [|I7|]-[[2"CJ|] show that 
the phase transition does not necessarily complete (as in old inflation) for some ranges of parame- 



ters of the GW stabilization mechanism. Ref. [17] found that the bubble nucleation rate r ~ e b is 
suppressed by the tunneling action 

A 2 M 2 e -2/v/£ -2/^3 

S p —= = (4.19) 



at temperature T, so the transition is slower if the radion mass nit is small or if e is large. Ref. [19] 
obtains quantitatively different but qualitatively similar results. 

An interesting recent development is the prediction of gravity waves produced during the phase 
transition, which can be large enough to be observable at LISA if the phase transition is strongly 
enough first order p0[]. The gravity waves are produced by bubble collisions and by turbulence. 
However they are only observably large for model parameters where perturbation theory is starting 
to break down, so it would be interesting to undertake a more careful calculation of the effect. 

5. The DGP model 

Dvali, Gabadadze and Porrati (DGP) found a remarkably simple way in which branes can 
modify gravity at large distances p(J], using the action 



S = 2M 2 p J d A x^f^R w + 2M| J d 4 xdy \/^G/? (5) (5. 1) 

where g and G is are the 4D and 5D metric determinants, respectively. By including the 4D Einstein 
term on a 3-brane, gravity is quasilocalized there, up to the distance scale 

Ml 
Ml 

At distances r < r c , gravity looks 4D, while for r > r c it looks 5D. The cosmology of this model 



was first studied in [27, 28], which found the modified Friedmann equation 

H 




(5.3) 



3Mj ' 4r ( 2 ^2r, ° A) 
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d. 



L 



{\ +z )fdx^- (5.6) 
Jo Hix) 



The upper sign leads to conventional cosmology as p — > 0, while the lower one has the interesting 
property that H — ► \/r c even as the matter density on the brane vanishes. This is called the self- 
accelerating solution, and gives rise to the possibility of explaining the accelerated expansion of 
the universe without an explicit cosmological constant term. To match the current Hubble rate, a 
small 5D gravity scale is required: 

M 5 = (M 2 p H Q ) 1/3 ~ 10 MeV (5.5) 

In principle this kind of acceleration is experimentally distinguishable from that due to dark 
energy, even if it has an arbitrary (but constant) equation of state w = p/p. The luminosity distance 
dL{z) as a function of redshift z, 

■H(x) 

has a different shape than in conventional dark energy models, so if the expansion history can be 
probed accurately as a function of z, as with type I supernovae, one can test the DGP model relative 
to fluid models. 

Interestingly, the DGP model also predicts nonstandard gravity at much shorter scales, poten- 
tially including the solar system. This comes about because of an extra scalar graviton polarization 
71 which changes the strength of the Newtonian gravitational potential in a distance-dependent 
way. One can think of % as a "brane-bending" mode, associated with transverse fluctuations of the 
brane. It modifies the strength of gravity at distances r > r* = (r^) 1 ' 3 , where r g = 2M/M 2 is the 
gravitational radius of a source of mass M, for example the sun if we are the testing solar system. 
Thus 

f2MMl\ 1/3 
* = {-Mf) 

Naively, the extra polarization % would spoil GR down to zero distance, but nonlinear terms in 
the Einstein equations become important when r < and these suppress the scalar contribution 
to the gravitational potential [g9|]. The different distance regimes and corresponding gravitational 
potentials are illustrated in fig. [2| 

GM 4GM _4GMr c 

r 3r 3r 2 



M r* r 



c 



Figure 2: Distance scales and corresponding behavior of gravitational potential in the DGP model. 

Of course even in the 4D regime, there are small corrections to the Newtonian potential, and 
these lead to anomalous precession of planetary orbits, and that of the moon. The precession rate 
turns out to be independent of the masses of the source or the orbiting body (see [30]), 

dkti 3 

-— = =F— = =F 5 ^ as/year (5.8) 
at 8ro 

which might be probed by future lunar laser ranging experiments. 
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There is yet another mass/distance scale associated with the % mode, 




A 



(10 3 



km) 



(5.9) 



For r < r\, the self coupling (dn) 2 \37i/ A 3 of the effective Lagrangian for the % becomes strong. 
If generic operators of the form (dn) 2N / A 4N ~ 4 were also present, as would be true for an ordinary 
Goldstone boson, this strong coupling would cause the theory to lose predictivity at any r < r\. 
Remarkably, such operators are not present in the DGP model, due to the unconventional shift 
symmetry d^n —> d^n + Cu, which is a remnant of the 5D Lorentz symmetry. Hence predictivity 
is not spoiled in the DGP model. However, precisely the absence of the (dn) 4 /A 4 term leads to a 
violation of causality and locality, as shown by ref. [f3l|]. This indicates that the DGP model has 
no ultraviolet completion which respects analyticity of the S -matrix, hence making it appear quite 
unlikely that it can be embedded within string theory. 
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